3.1.97 \(\int \frac {(d+e x^2)^2 (a+b \csc ^{-1}(c x))}{x^3} \, dx\) [97]

3.1.97.1 Optimal result
3.1.97.2 Mathematica [A] (verified)
3.1.97.3 Rubi [A] (verified)
3.1.97.4 Maple [A] (verified)
3.1.97.5 Fricas [F]
3.1.97.6 Sympy [F]
3.1.97.7 Maxima [F]
3.1.97.8 Giac [F]
3.1.97.9 Mupad [F(-1)]

3.1.97.1 Optimal result

Integrand size = 21, antiderivative size = 189 \[ \int \frac {\left (d+e x^2\right )^2 \left (a+b \csc ^{-1}(c x)\right )}{x^3} \, dx=-\frac {b c d^2 \sqrt {1-\frac {1}{c^2 x^2}}}{4 x}+\frac {b e^2 \sqrt {1-\frac {1}{c^2 x^2}} x}{2 c}+\frac {1}{4} b c^2 d^2 \csc ^{-1}(c x)+i b d e \csc ^{-1}(c x)^2-\frac {d^2 \left (a+b \csc ^{-1}(c x)\right )}{2 x^2}+\frac {1}{2} e^2 x^2 \left (a+b \csc ^{-1}(c x)\right )-2 b d e \csc ^{-1}(c x) \log \left (1-e^{2 i \csc ^{-1}(c x)}\right )+2 b d e \csc ^{-1}(c x) \log \left (\frac {1}{x}\right )-2 d e \left (a+b \csc ^{-1}(c x)\right ) \log \left (\frac {1}{x}\right )+i b d e \operatorname {PolyLog}\left (2,e^{2 i \csc ^{-1}(c x)}\right ) \]

output
1/4*b*c^2*d^2*arccsc(c*x)+I*b*d*e*arccsc(c*x)^2-1/2*d^2*(a+b*arccsc(c*x))/ 
x^2+1/2*e^2*x^2*(a+b*arccsc(c*x))-2*b*d*e*arccsc(c*x)*ln(1-(I/c/x+(1-1/c^2 
/x^2)^(1/2))^2)+2*b*d*e*arccsc(c*x)*ln(1/x)-2*d*e*(a+b*arccsc(c*x))*ln(1/x 
)+I*b*d*e*polylog(2,(I/c/x+(1-1/c^2/x^2)^(1/2))^2)-1/4*b*c*d^2*(1-1/c^2/x^ 
2)^(1/2)/x+1/2*b*e^2*x*(1-1/c^2/x^2)^(1/2)/c
 
3.1.97.2 Mathematica [A] (verified)

Time = 0.50 (sec) , antiderivative size = 194, normalized size of antiderivative = 1.03 \[ \int \frac {\left (d+e x^2\right )^2 \left (a+b \csc ^{-1}(c x)\right )}{x^3} \, dx=\frac {1}{4} \left (-\frac {2 a d^2}{x^2}+2 a e^2 x^2-\frac {2 b d^2 \csc ^{-1}(c x)}{x^2}+\frac {2 b e^2 x \left (\sqrt {1-\frac {1}{c^2 x^2}}+c x \csc ^{-1}(c x)\right )}{c}-\frac {b d^2 \left (-1+c^2 x^2+c^2 x^2 \sqrt {-1+c^2 x^2} \arctan \left (\sqrt {-1+c^2 x^2}\right )\right )}{c \sqrt {1-\frac {1}{c^2 x^2}} x^3}-8 b d e \csc ^{-1}(c x) \log \left (1-e^{2 i \csc ^{-1}(c x)}\right )+8 a d e \log (x)+4 i b d e \left (\csc ^{-1}(c x)^2+\operatorname {PolyLog}\left (2,e^{2 i \csc ^{-1}(c x)}\right )\right )\right ) \]

input
Integrate[((d + e*x^2)^2*(a + b*ArcCsc[c*x]))/x^3,x]
 
output
((-2*a*d^2)/x^2 + 2*a*e^2*x^2 - (2*b*d^2*ArcCsc[c*x])/x^2 + (2*b*e^2*x*(Sq 
rt[1 - 1/(c^2*x^2)] + c*x*ArcCsc[c*x]))/c - (b*d^2*(-1 + c^2*x^2 + c^2*x^2 
*Sqrt[-1 + c^2*x^2]*ArcTan[Sqrt[-1 + c^2*x^2]]))/(c*Sqrt[1 - 1/(c^2*x^2)]* 
x^3) - 8*b*d*e*ArcCsc[c*x]*Log[1 - E^((2*I)*ArcCsc[c*x])] + 8*a*d*e*Log[x] 
 + (4*I)*b*d*e*(ArcCsc[c*x]^2 + PolyLog[2, E^((2*I)*ArcCsc[c*x])]))/4
 
3.1.97.3 Rubi [A] (verified)

Time = 0.82 (sec) , antiderivative size = 228, normalized size of antiderivative = 1.21, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.238, Rules used = {5764, 5230, 27, 7293, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\left (d+e x^2\right )^2 \left (a+b \csc ^{-1}(c x)\right )}{x^3} \, dx\)

\(\Big \downarrow \) 5764

\(\displaystyle -\int \left (\frac {d}{x^2}+e\right )^2 x^3 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )d\frac {1}{x}\)

\(\Big \downarrow \) 5230

\(\displaystyle \frac {b \int -\frac {-\frac {d^2}{x^2}-4 e \log \left (\frac {1}{x}\right ) d+e^2 x^2}{2 \sqrt {1-\frac {1}{c^2 x^2}}}d\frac {1}{x}}{c}-\frac {d^2 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )}{2 x^2}-2 d e \log \left (\frac {1}{x}\right ) \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )+\frac {1}{2} e^2 x^2 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )\)

\(\Big \downarrow \) 27

\(\displaystyle -\frac {b \int \frac {-\frac {d^2}{x^2}-4 e \log \left (\frac {1}{x}\right ) d+e^2 x^2}{\sqrt {1-\frac {1}{c^2 x^2}}}d\frac {1}{x}}{2 c}-\frac {d^2 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )}{2 x^2}-2 d e \log \left (\frac {1}{x}\right ) \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )+\frac {1}{2} e^2 x^2 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )\)

\(\Big \downarrow \) 7293

\(\displaystyle -\frac {b \int \left (-\frac {d^2}{\sqrt {1-\frac {1}{c^2 x^2}} x^2}-\frac {4 e \log \left (\frac {1}{x}\right ) d}{\sqrt {1-\frac {1}{c^2 x^2}}}+\frac {e^2 x^2}{\sqrt {1-\frac {1}{c^2 x^2}}}\right )d\frac {1}{x}}{2 c}-\frac {d^2 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )}{2 x^2}-2 d e \log \left (\frac {1}{x}\right ) \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )+\frac {1}{2} e^2 x^2 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )\)

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {d^2 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )}{2 x^2}-2 d e \log \left (\frac {1}{x}\right ) \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )+\frac {1}{2} e^2 x^2 \left (a+b \arcsin \left (\frac {1}{c x}\right )\right )-\frac {b \left (-\frac {1}{2} c^3 d^2 \arcsin \left (\frac {1}{c x}\right )-2 i c d e \operatorname {PolyLog}\left (2,e^{2 i \arcsin \left (\frac {1}{c x}\right )}\right )-2 i c d e \arcsin \left (\frac {1}{c x}\right )^2+4 c d e \arcsin \left (\frac {1}{c x}\right ) \log \left (1-e^{2 i \arcsin \left (\frac {1}{c x}\right )}\right )-4 c d e \log \left (\frac {1}{x}\right ) \arcsin \left (\frac {1}{c x}\right )+\frac {c^2 d^2 \sqrt {1-\frac {1}{c^2 x^2}}}{2 x}-e^2 x \sqrt {1-\frac {1}{c^2 x^2}}\right )}{2 c}\)

input
Int[((d + e*x^2)^2*(a + b*ArcCsc[c*x]))/x^3,x]
 
output
-1/2*(d^2*(a + b*ArcSin[1/(c*x)]))/x^2 + (e^2*x^2*(a + b*ArcSin[1/(c*x)])) 
/2 - 2*d*e*(a + b*ArcSin[1/(c*x)])*Log[x^(-1)] - (b*((c^2*d^2*Sqrt[1 - 1/( 
c^2*x^2)])/(2*x) - e^2*Sqrt[1 - 1/(c^2*x^2)]*x - (c^3*d^2*ArcSin[1/(c*x)]) 
/2 - (2*I)*c*d*e*ArcSin[1/(c*x)]^2 + 4*c*d*e*ArcSin[1/(c*x)]*Log[1 - E^((2 
*I)*ArcSin[1/(c*x)])] - 4*c*d*e*ArcSin[1/(c*x)]*Log[x^(-1)] - (2*I)*c*d*e* 
PolyLog[2, E^((2*I)*ArcSin[1/(c*x)])]))/(2*c)
 

3.1.97.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 5230
Int[((a_.) + ArcSin[(c_.)*(x_)]*(b_.))*((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_ 
)^2)^(p_.), x_Symbol] :> With[{u = IntHide[(f*x)^m*(d + e*x^2)^p, x]}, Simp 
[(a + b*ArcSin[c*x])   u, x] - Simp[b*c   Int[SimplifyIntegrand[u/Sqrt[1 - 
c^2*x^2], x], x], x]] /; FreeQ[{a, b, c, d, e, f, m}, x] && NeQ[c^2*d + e, 
0] && IntegerQ[p] && (GtQ[p, 0] || (IGtQ[(m - 1)/2, 0] && LeQ[m + p, 0]))
 

rule 5764
Int[((a_.) + ArcCsc[(c_.)*(x_)]*(b_.))^(n_.)*(x_)^(m_.)*((d_.) + (e_.)*(x_) 
^2)^(p_.), x_Symbol] :> -Subst[Int[(e + d*x^2)^p*((a + b*ArcSin[x/c])^n/x^( 
m + 2*(p + 1))), x], x, 1/x] /; FreeQ[{a, b, c, d, e, n}, x] && IGtQ[n, 0] 
&& IntegerQ[m] && IntegerQ[p]
 

rule 7293
Int[u_, x_Symbol] :> With[{v = ExpandIntegrand[u, x]}, Int[v, x] /; SumQ[v] 
]
 
3.1.97.4 Maple [A] (verified)

Time = 4.46 (sec) , antiderivative size = 274, normalized size of antiderivative = 1.45

method result size
parts \(a \left (\frac {e^{2} x^{2}}{2}-\frac {d^{2}}{2 x^{2}}+2 d e \ln \left (x \right )\right )+i b d e \operatorname {arccsc}\left (c x \right )^{2}-\frac {b c \,d^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}{4 x}+\frac {b \,c^{2} d^{2} \operatorname {arccsc}\left (c x \right )}{4}-\frac {b \,\operatorname {arccsc}\left (c x \right ) d^{2}}{2 x^{2}}+\frac {b \,e^{2} \operatorname {arccsc}\left (c x \right ) x^{2}}{2}+\frac {b \,e^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, x}{2 c}-\frac {i b \,e^{2}}{2 c^{2}}-2 b d e \,\operatorname {arccsc}\left (c x \right ) \ln \left (1+\frac {i}{c x}+\sqrt {1-\frac {1}{c^{2} x^{2}}}\right )+2 i b d e \operatorname {polylog}\left (2, -\frac {i}{c x}-\sqrt {1-\frac {1}{c^{2} x^{2}}}\right )-2 b d e \,\operatorname {arccsc}\left (c x \right ) \ln \left (1-\frac {i}{c x}-\sqrt {1-\frac {1}{c^{2} x^{2}}}\right )+2 i b d e \operatorname {polylog}\left (2, \frac {i}{c x}+\sqrt {1-\frac {1}{c^{2} x^{2}}}\right )\) \(274\)
derivativedivides \(c^{2} \left (\frac {a \,x^{2} e^{2}}{2 c^{2}}+\frac {2 a d e \ln \left (c x \right )}{c^{2}}-\frac {a \,d^{2}}{2 c^{2} x^{2}}+\frac {i b d e \operatorname {arccsc}\left (c x \right )^{2}}{c^{2}}-\frac {b \,d^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}{4 c x}+\frac {b \,\operatorname {arccsc}\left (c x \right ) d^{2}}{4}-\frac {b \,\operatorname {arccsc}\left (c x \right ) d^{2}}{2 c^{2} x^{2}}+\frac {b \,\operatorname {arccsc}\left (c x \right ) x^{2} e^{2}}{2 c^{2}}+\frac {b \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, e^{2} x}{2 c^{3}}-\frac {i b \,e^{2}}{2 c^{4}}-\frac {2 b \ln \left (1+\frac {i}{c x}+\sqrt {1-\frac {1}{c^{2} x^{2}}}\right ) d e \,\operatorname {arccsc}\left (c x \right )}{c^{2}}-\frac {2 b \ln \left (1-\frac {i}{c x}-\sqrt {1-\frac {1}{c^{2} x^{2}}}\right ) d e \,\operatorname {arccsc}\left (c x \right )}{c^{2}}+\frac {2 i b \operatorname {polylog}\left (2, -\frac {i}{c x}-\sqrt {1-\frac {1}{c^{2} x^{2}}}\right ) d e}{c^{2}}+\frac {2 i b \operatorname {polylog}\left (2, \frac {i}{c x}+\sqrt {1-\frac {1}{c^{2} x^{2}}}\right ) d e}{c^{2}}\right )\) \(309\)
default \(c^{2} \left (\frac {a \,x^{2} e^{2}}{2 c^{2}}+\frac {2 a d e \ln \left (c x \right )}{c^{2}}-\frac {a \,d^{2}}{2 c^{2} x^{2}}+\frac {i b d e \operatorname {arccsc}\left (c x \right )^{2}}{c^{2}}-\frac {b \,d^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}{4 c x}+\frac {b \,\operatorname {arccsc}\left (c x \right ) d^{2}}{4}-\frac {b \,\operatorname {arccsc}\left (c x \right ) d^{2}}{2 c^{2} x^{2}}+\frac {b \,\operatorname {arccsc}\left (c x \right ) x^{2} e^{2}}{2 c^{2}}+\frac {b \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, e^{2} x}{2 c^{3}}-\frac {i b \,e^{2}}{2 c^{4}}-\frac {2 b \ln \left (1+\frac {i}{c x}+\sqrt {1-\frac {1}{c^{2} x^{2}}}\right ) d e \,\operatorname {arccsc}\left (c x \right )}{c^{2}}-\frac {2 b \ln \left (1-\frac {i}{c x}-\sqrt {1-\frac {1}{c^{2} x^{2}}}\right ) d e \,\operatorname {arccsc}\left (c x \right )}{c^{2}}+\frac {2 i b \operatorname {polylog}\left (2, -\frac {i}{c x}-\sqrt {1-\frac {1}{c^{2} x^{2}}}\right ) d e}{c^{2}}+\frac {2 i b \operatorname {polylog}\left (2, \frac {i}{c x}+\sqrt {1-\frac {1}{c^{2} x^{2}}}\right ) d e}{c^{2}}\right )\) \(309\)

input
int((e*x^2+d)^2*(a+b*arccsc(c*x))/x^3,x,method=_RETURNVERBOSE)
 
output
a*(1/2*e^2*x^2-1/2*d^2/x^2+2*d*e*ln(x))+I*b*d*e*arccsc(c*x)^2-1/4*b*c*d^2/ 
x*((c^2*x^2-1)/c^2/x^2)^(1/2)+1/4*b*c^2*d^2*arccsc(c*x)-1/2*b*arccsc(c*x)* 
d^2/x^2+1/2*b*e^2*arccsc(c*x)*x^2+1/2*b/c*e^2*((c^2*x^2-1)/c^2/x^2)^(1/2)* 
x-1/2*I*b/c^2*e^2-2*b*d*e*arccsc(c*x)*ln(1+I/c/x+(1-1/c^2/x^2)^(1/2))+2*I* 
b*d*e*polylog(2,-I/c/x-(1-1/c^2/x^2)^(1/2))-2*b*d*e*arccsc(c*x)*ln(1-I/c/x 
-(1-1/c^2/x^2)^(1/2))+2*I*b*d*e*polylog(2,I/c/x+(1-1/c^2/x^2)^(1/2))
 
3.1.97.5 Fricas [F]

\[ \int \frac {\left (d+e x^2\right )^2 \left (a+b \csc ^{-1}(c x)\right )}{x^3} \, dx=\int { \frac {{\left (e x^{2} + d\right )}^{2} {\left (b \operatorname {arccsc}\left (c x\right ) + a\right )}}{x^{3}} \,d x } \]

input
integrate((e*x^2+d)^2*(a+b*arccsc(c*x))/x^3,x, algorithm="fricas")
 
output
integral((a*e^2*x^4 + 2*a*d*e*x^2 + a*d^2 + (b*e^2*x^4 + 2*b*d*e*x^2 + b*d 
^2)*arccsc(c*x))/x^3, x)
 
3.1.97.6 Sympy [F]

\[ \int \frac {\left (d+e x^2\right )^2 \left (a+b \csc ^{-1}(c x)\right )}{x^3} \, dx=\int \frac {\left (a + b \operatorname {acsc}{\left (c x \right )}\right ) \left (d + e x^{2}\right )^{2}}{x^{3}}\, dx \]

input
integrate((e*x**2+d)**2*(a+b*acsc(c*x))/x**3,x)
 
output
Integral((a + b*acsc(c*x))*(d + e*x**2)**2/x**3, x)
 
3.1.97.7 Maxima [F]

\[ \int \frac {\left (d+e x^2\right )^2 \left (a+b \csc ^{-1}(c x)\right )}{x^3} \, dx=\int { \frac {{\left (e x^{2} + d\right )}^{2} {\left (b \operatorname {arccsc}\left (c x\right ) + a\right )}}{x^{3}} \,d x } \]

input
integrate((e*x^2+d)^2*(a+b*arccsc(c*x))/x^3,x, algorithm="maxima")
 
output
1/2*a*e^2*x^2 + 1/4*b*d^2*((c^4*x*sqrt(-1/(c^2*x^2) + 1)/(c^2*x^2*(1/(c^2* 
x^2) - 1) - 1) - c^3*arctan(c*x*sqrt(-1/(c^2*x^2) + 1)))/c - 2*arccsc(c*x) 
/x^2) + 2*a*d*e*log(x) - 1/2*a*d^2/x^2 + 1/4*(4*I*b*c^2*d*e*log(-c*x + 1)* 
log(x) + 4*I*b*c^2*d*e*log(x)^2 + 4*I*b*c^2*d*e*dilog(c*x) + 4*I*b*c^2*d*e 
*dilog(-c*x) + 2*(b*c^2*arctan2(1, sqrt(c*x + 1)*sqrt(c*x - 1)) + I*b*c^2* 
log(c))*e^2*x^2 + I*b*e^2*log(c*x - 1) - I*(b*e^2*(log(c*x + 1)/c^2 + log( 
c*x - 1)/c^2) + 16*b*d*e*integrate(1/2*log(x)/(c^2*x^3 - x), x))*c^2 + 4*c 
^2*integrate(1/2*(b*e^2*x^2 + 4*b*d*e*log(x))*sqrt(c*x + 1)*sqrt(c*x - 1)/ 
(c^2*x^3 - x), x) + (-I*b*c^2*e^2*x^2 - 4*I*b*c^2*d*e*log(x))*log(c^2*x^2) 
 + (4*I*b*c^2*d*e*log(x) + I*b*e^2)*log(c*x + 1) - 2*(-I*b*c^2*e^2*x^2 - 4 
*(b*c^2*arctan2(1, sqrt(c*x + 1)*sqrt(c*x - 1)) + I*b*c^2*log(c))*d*e)*log 
(x))/c^2
 
3.1.97.8 Giac [F]

\[ \int \frac {\left (d+e x^2\right )^2 \left (a+b \csc ^{-1}(c x)\right )}{x^3} \, dx=\int { \frac {{\left (e x^{2} + d\right )}^{2} {\left (b \operatorname {arccsc}\left (c x\right ) + a\right )}}{x^{3}} \,d x } \]

input
integrate((e*x^2+d)^2*(a+b*arccsc(c*x))/x^3,x, algorithm="giac")
 
output
integrate((e*x^2 + d)^2*(b*arccsc(c*x) + a)/x^3, x)
 
3.1.97.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\left (d+e x^2\right )^2 \left (a+b \csc ^{-1}(c x)\right )}{x^3} \, dx=\int \frac {{\left (e\,x^2+d\right )}^2\,\left (a+b\,\mathrm {asin}\left (\frac {1}{c\,x}\right )\right )}{x^3} \,d x \]

input
int(((d + e*x^2)^2*(a + b*asin(1/(c*x))))/x^3,x)
 
output
int(((d + e*x^2)^2*(a + b*asin(1/(c*x))))/x^3, x)